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Abstract 

This paper is devoted to estimates of the exponential decay of eigen- 
functions of difference operators on the lattice Z™ which are discrete 
analogs of the Schrodinger, Dirac and square-root Klein-Gordon opera- 
tors. Our investigation of the essential spectra and the exponential decay 
of eigenfunctions of the discrete spectra is based on the calculus of the 
so-called pseudodifference operators (i.e., pseudodifferential operators on 
the group Z n with analytic symbols, as developed in the paper [29]), and 
the limit operators method (see |34] and the references cited there). We 
obtain a description of the location of the essential spectra and estimates 
of the eigenfunctions of the discrete spectra of the main lattice opera- 
tors of quantum mechanics, namely: matrix Schrodinger operators on Z n , 
Dirac operators on Z 3 , and square root Klein-Gordon operators on Z n . 
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1 Introduction 

The main goal of this paper is estimates of the exponential decay of eigenfunc- 
tions of difference operators on the lattice Z™. Particular attention is paid to 
the discrete analogs of the Schrodinger, Dirac and square-root Klein-Gordon op- 
erators. Schrodinger operators on the lattice Z™ appear at several places, so in 
the tight binding model in solid state physics, in the propagation of spin waves 
and waves in quasi-crystals, and in mathematical models of nano-structure, to 
mention only a few (see, [3TJ [551 133 SI] ) • There is an extensive bibliography 
devoted to different aspects of the spectral theory of discrete Schrodinger op- 
erators. Let us only mention the papers [3J [H [T31 [TH [331 EH Hi], and see also 
the references cited there. In our recent papers [301 131] > we study the essential 
spectra of discrete Schrodinger operators with variable magnetic and electric 
potentials on the lattice Z ra and on periodic combinatorial graphs by means of 
the limit operators method (for the latter see also [51]). 
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In the last time, also discrete relativistic operators attracted many attention. 
They were used, e.g., in comparative studies of relativistic and nonrelativistic 
electronlocalization phenomena [5j , in relativistic investigations of electrical con- 
duction in disordered systems [35] i m the construction of supertransparent mod- 
els with supersymmetric structures |41j . and in relativistic tunnelling problems 

The problem of exponential estimates of solutions of the elliptic partial differ- 
ential equations with applications to the Schrodinger operator is a classical one. 
There are an extensive bibliography devoted to this problem (see for instance 
[IJ El El E4 EH EI])- Exponential estimates of solutions of pseudodifferential 
equations are considered in EHl Qj3 EH1 1231 ESI EI] • Note also the our recent 
papers [32l [33] where we proposed a new approach to exponential estimates for 
partial differential and pseudodifferential operators, based on the limit operators 
method (see [34]). 

Our approach to study essential spectra and the exponential decay of eigen- 
functions is based on the calculus of pseudodifference operators (i.e., pseudodif- 
ferential operators on the group Z") with analytic symbols as developed in [25] , 
and on the limit operators method (see [34] and the references cited there). 

The paper is organized as follows. In Section 2 we recall some auxiliary facts 
on the pseudodifference operators with analytic symbols on Z n , limit operators, 
essential spectra and the behavior of solutions of pseudodifference equations at 
infinity. 

In Section 3 we consider the discrete Schrodinger operators on Z 2 (Z™, C N ) 
of the form 

n 

(Hu)(x) = (V ek - e M ^>) (V_ efc - e -**<»)) u (x) + $(x)u(x) 
fc=i 

where V eh is the operator of shift by , the are real- valued bounded slowly 
oscillating functions on Z™ , and $ is a Hermitian slowly oscillating and bounded 
matrix function on Z™. We show that the essential spectrum sp ess H of H is 
the interval 

n 

sp ess H=\J[Xf,X^ + 4n] 

3 = 1 

where 

Xf l := liminf Aj($(a;)), Xf p := limsup Aj(*(a;)) 

x >oo x — > 00 

and where the Xj(^(x)) are the increasingly ordered eigenvalues of the matrix 
$(x), i.e. 

Ai($(a:)) < A 2 ($(x)) < . . . < X N (M.x)) 

for x £ Z" large enough. Note that sp ess H does not depend on the exponents 
afe, and that there is a gap (A^ up + 4n, X™^) in the essential spectrum of H if 
r ,; ' • i„ • A;; ; . 

We also obtain the following estimates of eigenfunctions belonging to points 
in the discrete spectrum of H . In each of the cases 
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A € (X" up + 4n, Xj+i) is an eigenvalue of H and 



< r < cosh 



_j , mm{A - A™ p - 2rc, Aj^ - A + 2n}" 



2n 



• A > A^ lp + 4n is an eigenvalue of H and 

A - A s " p - 2n 



< r < cosh 



-1 I n 



2n 



A < A5" f is an eigenvalue of H and 

_! /Ai nf -A + 2n 



< r < cosh 



2n 



every A-eigenfunction uoiH has the property that e^'it € l p (Z n , C N ) for every 
1 < p < oo. 

In Section 4 we introduce self-adjoint Dirac operators on the lattice 1? with 
variable slowly oscillating electric potentials. In accordance with the general 
properties of Dirac operators on K 3 (see for instance 0142]), the corresponding 
discrete Dirac operator on Z 3 should be a self-adjoint system of the first order 
difference operators. We thus let 

V := V + e<PE 4 (1) 

where 

£>o := chd k -f k + c 2 m-f°, 

E N is the N x N unit matrix, the 7 fe with k = 0, 1, 2, 3 refer to the 4x4 Dirac 
matrices, the 

d k :=I-V eh , k= 1,2,3 

are difference operators of the first order, h is Planck's constant, c is the speed of 
light, m and e are the mass and the charge of the electron and, finally, $ is the 
real electric potential. The operator T>, acting on Z 2 (Z 3 , C 4 ), can be considered 
as the direct discrete analog of the Dirac operator on K 3 , but note that T> is 
not self-adjoint on l 2 (Z, 3 , C 4 ). To force the self-adjointness, we consider the 
"symmetrization" B := D + e$J of V with 



V 
VI 

which acts on l 2 (Z 3 , C 8 ). The operator D is self-adjoint, and 
n2 _( (h 2 c 2 T + m 2 c*)E 4 



(2) 



(h 2 c 2 T + m 2 c 4 )Ei 
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where h 2 c 2 T + m 2 c 4 is the lattice Klein-Gordon Hamiltonian with Laplacian 

3 3 

T:=J2d* k d k = J2(^-V ek -V: k ). 



fc=i k=i 



Note that one-dimensional Dirac operators of the form @ on Z were considered 

in gaisg. 

We prove that the essential spectrum of D is the union 



sPess D = [ e$hlf - \/l2ft 2 c 2 + m 2 c 4 , e$ sup - mc 2 ] 



U [e$ inf + mc 2 , e$ sup + v / 12fi 2 c 2 + m 2 c 4 ], 

where 



<3> := Uminf $ sup := limsup$(x). 



.- lllllllll Y - .= J 

a; >oo ^ — > qq 

Again we observe that if e<I> sup — e<I> mf < 2mc 2 , then the essential spectrum of 
M has the gap (e$ sup - mc 2 , e$ inf + mc 2 ). 

We also obtain the following estimates of eigenfunctions of the discrete spec- 
trum. Let A be a point of the discrete spectrum, and let A and r > satisfy one 
of the conditions 

• A e (e$ sup - mc 2 , e<I> inf + mc 2 ) and 



< r < cosh 1 



{(e$ inf - A) 2 , (e$ sup - A) 2 } + 6ft 2 c 2 
6 ft 2 c 2 



• A > e$ sup + Vl2ft 2 c 2 + m 2 c 4 and 

_ x / (e$ sup - A) 2 - m 2 c 4 - 6ft 2 c 2 



< r < cosh 



6ft 2 c 2 



A < e$ inf - v / 12ft 2 c 2 + m 2 c 4 and 

, / (e$ inf - A) 2 - m 2 c 4 - 6ft 2 c 2 
< r < cosh ' 



Then every A-eigenfunction u of the operator D satisfies e r ^u £ l p (l?, C 8 ) for 
every p 6 (1, oo). 

In Section 5, we consider the lattice model of the relativistic square root 
Klein-Gordon operator as the pseudodiffcrencc operator of the form 



K. := \Jc 2 h 2 Y + m 2 c 4 + e$ 

on l 2 (X n ). We determine the essential spectrum of JC and obtain exact estimates 
of the exponential decay at infinity of eigenfunctions of the discrete spectrum. 
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2 Pseudodifference operators, essential spectra, 
and exponential estimates 

2.1 Some function spaces 

For each Banach space X, B(X) refers to the Banach algebra of all bounded 
linear operators acting on X. For 1 < p < oo, we let Z P (Z", C^) denote the 
Banach space of all functions on Z™ with values in with the norm 

ll/llf»(z»,c~) := J2 \\f( x )\\ P C N < oo if p < oo, 
||/||j°°(z«,c") := sup ||/(a;)|| C iv < oo. 

rrGZ™ 

The choice of the norm on C N is not of importance in general; only for p = 2 
we choose the Euclidean norm (such that ? 2 (Z™, C^) becomes a Hilbert space 
and B(C N ) a C*-algcbra in the usual way). Given a positive function w on Z™, 
which we will call a weig ht, let P(Z n , C N , w) stand for the Banach space of all 
functions on Z™ with values in such that 

IMI(p(z>\C N ,to) := ||^w|bp(z«,C JV ) < oo- 

Similarly, we write l°°(Z n , B(C N )) for the Banach algebra of all bounded func- 
tions on Z n with values in B(C N ) and the norm 

ll/lh~(z«,s(c N )) : = SU P ll/(z)lls(C w ) <°°- 

xSZ" 

Finally, we call a function a G ?°°(Z ra , ^(C^)) slowly oscillating if 

lim \\a(x + y)- a(a;)||B(c«) = 

for every point y £ Z" . We denote the class of all slowly oscillating functions 
by SO(Z n , B{C N )) and write simply SO(Z n ) in case N = 1. 

2.2 Pseudodifference operators 

Consider the n-dimensional torus T" as a multiplicative group and let 

_ / 1 \ n (iti ■...■*„ _ / 1 
M ' V2^y ti-...-t n ~ yhd) T 

denote the corresponding normalized Haar measure on T™. 

Definition 1 Let S(N) denote the class of all matrix-valued functions a = 
( a ij)i,j=i on z ™ x T ™ with 

\\a\\ k := sup ||a t Q a(x,t)|| B(c « ) < oo (3) 

(l,t)£Z"xT», \a\<k 
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for every non-negative integer k, provided with the convergence defined by the 
semi-norms \a\k- To each function a € S(N), we associate the pseudodifference 
operator 

(Op(o)«)(x) := f a(x,t)u(t)t x dfi(t), xeZ n , (4) 

which is defined on vector-valued functions with finite support. Here, u refers 
to the discrete Fourier transform of u, i.e., 

u(t) := J2 u{x)t x , t e T\ 

We denote the class of all pseudodifference operators by OPS(N). 

Pseudodifference operators on Z™ can be thought of as the discrete analog 
of pseudodifferential operators on K ra (see for instance [40l 03]); they can be 
also interpreted as (abstract) pseudodifferential operators with respect to the 
group Z™. For another representation of pseudodifference operators, we need 
the operator V a of shift by a E Z n , i.e. the operator V a on l p (Z n , C ) which 
acts via 

(V a u)(x) = u(x -a), x £ Z n . 
Then the operator Op (a) can be written as 

Op (a) = a ^ ( 5 ) 

where 

a a (x) := I a(x,t)t a dfi(t). 
Integrating by parts we obtain 

l|aa||i»(Z»,B(c")) < C|a|2(l + \u\)~ 2 , (6) 

whence 

l|Op(a)|| w(z „ iC jv) := ^2 ll a alli°°(Z",B(C«)) < 00 ■ ( 7 ) 

QgZ" 

We thus obtain that the pseudodifference operator Op (a) belongs to the Wiener 
algebra W(Z n ,C ) which, by definition, consists of all operators of the form 
(|5|) with norm |(7J). It is an immediate consequence of this fact that all operators 
Op (a) in OPS(N) are bounded on P(Z n ,C N ) for all p e [l,oo]. Moreover, 
since the algebra W(Z n ,C ) is inverse closed in B(l p CZ n ,C N )), the spectrum 
of Op (a) € OPS(N) is independent of the underlying space l p (Z n ,C N ). 
Note that the operator (jU can be also written as 

Op(o)u(x)= V / a(x,t)t x - y u{y)d f i(t), 
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which leads to the following generalization of pseudodiffcrence operators. Let 
a be a function on Z n x 2" x T" with values in 0(C) which is subject to the 
estimates 

\a\k=- sup \\d?a(x,y,t)\\ B(C N } < oo (8) 

(x,i/,t)eZ"xZ"xr ; | ct| < fc 

for every non- negative integer k. Let Sd(N) denote the set of all functions with 
these properties. To each function a 6 Sd(N), we associate the pseudodifference 
operator with double symbol 

(Op d (a)u)(x):= V f a(x,y,t)u(y)t x -y dfi(t) (9) 

where u : Z™ — > C w is a function with finite support. The right-hand side of ([9]) 
has to be understood as in (5.6) in [34], which is in analogy with the definition 
of an oscillatory integral (see [3D] and also Section 4.1.2 in [33] )• The class of 
all operators of this form is denoted by OPSd{N). 

The representation of operators on Z™ as pseudodifference operators is very 
convenient due to the fact that one has explicit formulas for products and ad- 
joints of such operators. The basic results are as follows. 

Proposition 2 (i) Let a, b £ S(N). Then the product Op(a)Op(&) is an op- 
erator in OPS(N), and Op (a) Op (b) = Op (c) with 

c(x,t)=y2 [ a(x,tr)b{x + y,T)T-y d/*(r), (10) 

with the right-hand side understood as an oscillatory integral. 

(ii) Let a £ S(N) and consider Op (a) as acting on l p (Z n , C N ) with p £ (1, oo). 
Then the adjoint operator of Op (a) belongs to OPS(N), too, and it is of the 
form Op (a)* = Op (6) with 

b(x,t)= V/ a*(x + y,tT)r- y d(i(r), (11) 

where a*{x,t) is the usual adjoint (i.e., transposed and complex conjugated) 
matrix. 

(Hi) Let a £ S d (N). Then Op d (a) £ OPS(N), and Op d (a) = Op (a*) where 
aF(x,£)= 2_, I a(x + y,tr)T~ v dfi(r). 

2.3 Limit operators and the essential spectrum 

Recall that an operator A £ B(X) is a Fredholm operator if its kernel kerA = 
{x £ X : Ax = 0} and its cokernel cokerA = X/(AX) are finite-dimensional 
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linear spaces. The essential spectrum of A consists of all points A e C such that 
the operator A — XI is not a Fredholm operator. We denote the (usual) spectrum 
and the essential spectrum of A by spec x ^4 and sp ess xA, respectively. 

Our main tool to study the Fredholm property is limit operators. The fol- 
lowing definition is crucial in what follows. 

Definition 3 Let A e B{P(Z n , C N )) with p <= (1, oo), and let h : N -> Z™ be a 
sequence which tends to infinity in the sense that \h(n)\ — ■> oo as n —* oo. An 
operator A h € B(l p (Z n , C N j) is called a limit operator of A with respect to the 
sequence h if 

s-lim m ^ 00 V r _ /l(m) AV r /l(m) = A h and s-lim m ^ 00 V_ /l(m )A*V r /l(m ) = (A h )* , 

where s-lim refers to the strong limit. Clearly, every operator has at most one 
limit operator with respect to a given sequence. We denote the set of all limit 
operators of A by op(A). 

Let al be the operator of multiplication by the function a € Z°°(Z", B(C N )). 
A standard Cantor diagonal argument shows that every sequence h tending to 
infinity possesses a subsequence g such that, for every ieZ", the limit 

lim a(x + g{m)) =: a 9 {x) 

m — >oo 

exists. Clearly, a 9 is again in l°°(7, n ,B(C N )). Hence, all limit operators of al 
are of the form a 9 1. In particular, if a € SO(1i n ,B(C N )), then it follows easily 
from the definition of a slowly oscillating function that all limit operators of al 
are of the form a 9 1 where now a 9 e B(C N ) is a constant function. 

Let Op (a) <G OPS(N), and let h : N — > Z" be a sequence tending to infinity. 
Then V_ h ( m )AV h ( m ) — Op (o m ) with a m {x) := a(x+h(m), t). It follows as above 
that the sequence h has a subsequence g such that a(x + g(m), t) converges to a 
limit a 9 (x, t) for every ieZ" uniformly with respect to t £ T". One can prove 
that the so-defined function a 9 belongs to S(N) and the associated operator 
Op (a 9 ) is the limit operator of Op (a) with respect to g. 

The following theorem gives a complete description of the essential spectrum 
of pseudodifference operators in terms of their limit operators. 

Theorem 4 Let a <G S(N). Then, for every p G (1, oo), 

s PessiP Op(a) = |J spec ; ,Op(a 9 ) (12) 

Op {as)£op(A) 

where re [l,oo] is arbitrary. 

Since spec ir Op(a s ) does not depend on the underlying space, the essential 
spectrum sp ess ; P Op (a) is independent of p e (1, oo). Hence, in what follows we 
will omit the explicit notation of the underlying space in the spectrum and the 
essential spectrum. 
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2.4 Pseudodifference operators with analytic symbols and 
exponential estimates of eigenfunctions 

For r > 1 let K r be the annulus {t G C : < |i| < r}, and let K™ be the 
product K r x . . . x K r of n factors. 

Definition 5 Let S(N, IK") denote the set of all functions 

a : Z n x K™ -» B(C*) 
which are analytic with respect to t in the domain K™ and satisfy the estimates 
|a| fe := V sup ||a t a a(a;,t)||B (C iv ) < oo 

a|</c r 

for every non-negative integer k. With every function a G S(N, K"), we asso- 
ciate a pseudodifference operator defined on vector-valued functions with finite 
support via Q , and we denote the corresponding class of pseudodifference op- 
erators by OPS{N, K"). 

Definition 6 For r > 1, let W(K") denote the class of all exponential weights 
w = expv, where v is the restriction onto Z" of a function v S C^(R n ) with 
the property that, for every point x G K" and every j = 1, n, 

dv( x i 

-logr < — ^ < logr. (13) 

dXj 

In what follows we will denote both the function v on K™ and its restriction 
onto Z™ by v. Note that it is an immediate consequence of Definition [5] that if 
w € W(K"), then w» G W(K£) for every /x G [-1, 1]. 

Proposition 7 Let A := Op (a) G OPS (AT, K") and w G W{K r r l ). Then the 
operator A w := wAw -1 , defined on vector-valued functions with finite support, 
belongs to the class OPSd(N), and A w = Op d (6) with 

b(x,y,t) =a{x,e~ e ™ (x < v) ■ t) 

where 

e -d w (x, y ) . t ._ L-e w ,i(^y) tlj e - 9 ™' 2( - x ' y h 2 , • ■ • , e- ew < n{x ' v) t^\ 

and 

f 1 dv(Sl-j)x+jy) 
°w,j(x, y) ^ try. 

Proposition [5] and estimate (| 1 5|) imply the following theorem. 

Theorem 8 Let a G S(N,K^) and w G W(K"). Then Op (a) is a bounded 
operator on each of the spaces Z P (Z™, C^, w) with 1 < p < 00. 
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Next we consider essential spectra of pseudodiffcrence operators on weighted 
spaces. Let a, A and A w be as in Proposition [7] One can easily check that for 
h G Z n 

V- h A w V h = Op d (b h ) with b h (x,y,t) = a(x + M^"^"^" 1 • t)). 

Let now /i : N -> Z™ be a sequence tending to infinity. Then there exists a 
subsequence g of h such that the limit operator of A w with respect to g exists 
and 

A 9 w = Op d (b 9 ) with b 9 (x,y,t) = a 9 {x,e e ^ x ' y) ■ t) (14) 

where 

a 9 (x,t):— lim a(x + g(m),t) (15) 

m— >oo 

and ] 

0&(a,y):= lim / Vw((l-7)a; + 7t/ + g(TO))d7. (16) 

The limits in (|15p and (fT!)]) are understood as pointwise with respect to x, y G Z" 
and uniform with respect to i G T™. 

Theorem 9 Let a G S(N,K?) and w G W(Kp) 7 set A := Op (a) and A w := 
iDylw" 1 , and consider A as operating from V (Z™ , C N , w) to Z P (Z™, C , w) where 
p G (1, oo). Then 

sp ess Op (a) = [J spec Op d (b 9 ) 

o Pd (6s)eop(A„) 

with b 9 as in \14\ - 

Remark 10 Note that the essential spectrum of an operator in OPS(N,¥L™), 
considered as acting on P(Z", C N , w), is independent of p G (l,oo), but it can 
depend on the weight in general. But if the weight w — e v has the property that 

lim Vv(x) = 0, (17) 

X — >oc 

then the symbol a 9 w does not depend on the weight and, hence, the essential 
spectrum of Op (a) G B(l p (1/ 1 ,C N ,w)) with a G 5(jV,K™) is independent both 
on p G (1, oo) and on the weight w. Important examples of weights satisfying 
( |_?7p are the power weights w{x) = (1 + |a;| 2 ) S//2 = ei log ( 1+ \ x \ ) with s > and 
the subexponential weights w(x) = e"' 1 ^ where a > and (3 G (0, 1). 

The next theorem provides exponential estimates of solutions of pseudodif- 
ference equations. 

Theorem 11 Let A = Op(a) G OPS(N,K r >) and w G W(K"). Suppose that 
lim-r^oo w(x) = +oo and that is not in the essential spectrum of A w ^ : 
P(Z n ,C N ) -► P(Z n ,C N ) for somep G (l,oo) and every fj, G [-1,1]. If u € 
l p (Z n , C , w _1 ) is a solution of the equation Au — f with f G / P (Z™, C , w), 
thenu(El p (Z n ,C N ,w). 
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Theorem 1111 has some important corollaries. 

Theorem 12 Let a, A and w be as in the previous theorem, and let A be an 
eigenvalue of A which is not in the essential spectrum of A w ^ : l p (Z n ,C ) — * 
l p (Z n , C N ) for somep £ (1, oo) and every \x £ [0, 1]. Then every X-eigenf unction 
belongs to l p (Z n , C , w) for every p £ (1, oo). 

Corollary 13 Let A — Op(a) £ OPS(N, K") and let A be an eigenvalue of 
A which is not in the essential spectrum of A. Then every X-eigenf unction 
u = (ui, ...,un) satisfies the sub- exponential estimate 

sup\ Ul (x)\ <C ie ~ alx ^ , x£Z n ,i = l,...,N (18) 
for arbitrary a > and < (3 < 1. 

Proof. Let w(x) = e v ( x > where v(x) — a\x\@ with a > and < < 1. Then 
liirix^oo Vn(i) = 0, whence A 9 wii = A 9 for every limit operator A 9 . Let A be an 
eigenvalue of A which is not in the essential spectrum of A. Then A is not in 
the essential spectrum of A w ^ for every /x £ [0, 1]. Hence, by Theorem [T2l every 
A-eigenfunction belongs to each of the spaces Z P (Z™,C , w) with p £ (l,oo). 
Applying the Holder inequality we obtain estimate (fT8| . ■ 

We are now going to specialize these results to the context of slowly oscil- 
lating symbols and slowly oscillating weights. 

Definition 14 The symbol a £ S(N, K") is said to be slowly oscillating if 

lim sup \\a(x + y,t) — a(x, t)||e(c«) = 
x->oo teT „ 

for every y £ Z™. We write S sl (N,K™) for the class of all slowly oscillating 
symbols and OPS sl (N,K.^:) for the corresponding class of pseudodifference op- 
erators. 

Definition 15 The weight w = e v £ W(K r ) is slowly oscillating if the partial 
derivatives are slowly oscillating for j = 1, . . . , n. We denote the class of 

all slowly oscillating weights by W sl (K r ). 

Example 16 Ifv(x) = j\x\, then — 7§r for j = 1, . . . , n. Thus, w := e" 
is in W sl (K r ) ifj<r. 

The next theorem describes the structure of the limit operators of the oper- 
ator A w = wAw- 1 if A e OPS sl {N,K?) and w £ W sl {K r ). 

Theorem 17 Let A = Op (a) £ OPS sl (N,K?) and w £ W sl (K r ). Then the 
limit operator A 9 ^ of A w with respect to the sequence g tending to infinity exists 
if the limits 

a g (t)= lim a(s(m),*), 6 9 w = lim (Vv)(g(m)) (19) 

m — >oo m — >oc 

exist. In this case, it is of the form 

A 9 W = Op (c g ) with c g {x,t) = a g {6 9 w -t). (20) 
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Consequently, if A and w are as in this theorem, then the limit operators A 9 W » 
are invariant with respect to shifts. This fact implies the following explicit de- 
scription of their essential spectra. Let {Aj(j4^ M )(f)}™ =1 denote the eigenvalues 
of the matrix a g {6 9 u ^ ■ t). Then 

N 

spec /P(z „ jC iY)A^ = \J{Xj{A 9 wfi )(t) :teT n and j = l,...,n}, 

3=1 

whence 

iV 

si?essir(zn,c»)Aw» = \J \J {\j(A 9 w „)(t) : t e T 1 and j = 1, n} . 

3 Matrix Schrodinger operators 
3.1 Essential spectrum 

In this section we consider the essential spectrum and the behavior at infin- 
ity of eigenfunctions of general discrete Schrodinger operators acting on u e 
l 2 (Z n ,C N ) by 

n 

(Hu)(x) = ^(V efe - e <0 * (x) ) (V- e „ - e- ia «W)u(x) + $(x) (21) 
fe=i 

where = (0, . . . , 0, 1, 0, . . . , 0) with the 1 standing at the fcth place, the 
a k e SO(Z n ) are real-valued, and $ e SO(Z n ,B(C N )) is Hcrmitian. The vec- 
tor a := (oi, . . . , a n ) is the discrete analog of the magnetic potential, whereas 
$ can be viewed of as a discrete analog of the electric potential. Since the es- 
sential spectrum of H is independent of p G (1, oo), we consider the case p = 2 
only. Note that our assumptions guarantee that H is a self-adjoint operator on 
l 2 {Z n ,C N ). 

The limit operators H g of H are of form 

n 

H 9 = 52{V ek -e<I){V- ek -e-<l) + &I 

k=l 

n 

= £)(2J-e- io ^_ et -e ia fy ek ) + &I 
k=i 

with the constant functions 

a£ = lim a,k{x + g{m)) and $ s = lim $>(x + g(m)). 

m — >oc m — >oo 

Let U : l 2 (Z n ,C N ) -» ^(Z™,^) be the unitary operator 

(t/u)(z) = e-'^^uCx), a 9 = (a 9 ,...,a 9 n ). 
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Then 

n 

U*H g U = Y,(2I-V- ek -V eh )+&. 

k=l 

Further, the operator H' g := U*H g U is unitarily equivalent to the operator of 
multiplication by the function 



fc=i 

acting on L 2 {[0, 2<ir] n , C N ). Hence, 



Hgtyx, . . . , Vn) := 4 £ sin 2 ^ + $9, V'fe € [0, 2tt], 
fc=i 



JV 

spec Bp = specif = |J [Aj($ 9 ), A,($ ff ) + 4n], 

3=1 

where the Aj($ 9 ) refer to the eigenvalues of the matrix $ 9 . Applying formula 
(JT3J) we obtain 

JV 

sPess H = \J{J [A,-(<F), Xji® 9 ) + in] (22) 
g 3=1 

where the first union is taken over all sequences g for which the limit operator 
of H exists. Let Aj($(x)), j — 1, . . . , N, denote the eigenvalues of the matrix 
We suppose that these eigenvalues are simple for x large enough and that 
they are increasingly ordered, 

Ai($(a:)) < Aa($(a;)) < . . • < A* (*(«)). 

Then one can show that the functions x i— > Aj($(a;)) belong to SO(Z n ). Let 

Af f := liminf Aj($(a;)), A™ p := limsup A,-(#(a:)). 

x >oo x — >DO 

Since the set of the partial limits of a slowly oscillating function on Z™ is con- 
nected for n > 1 (see [M], Theorem 2.4.7), we conclude from that 

JV 

sp ess i?=|J[Af ) Af p + 4n] (23) 
j=i 

for n > 1 . Note that if A* up + 4n< Aj 1 ^ , then there is the gap (A" up + An, Aj^ f x ) 
in the essential spectrum of H . 

In case n = 1, the set of the partial limits of a slowly oscillating function 
on Z consists of two connected components, which collect the partial limits as 
x — > — oo and x — > +oo, respectively. Accordingly, in this case we set 

Af f,± := Hminf \j{${x)), A* up,± := limsup A,-(* (a)). 

and obtain 

>!>.,, // (J ([Af '-, Af p < + 4] U [Af •+ Af p '+ + 4] 

3 = 1 
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3.2 Exponential estimates of eigenfunctions 

Our next goal is to apply Theorem [TT] to eigenfunctions of (discrete) eigenval- 
ues of the operator H with slowly oscillating potentials. We will formulate the 
results for n > 1 only; for n = 1 the non-connectedness of the set of the par- 
tial limits requires some evident modifications. According to (j2li|) , the discrete 
spectrum of H is located outside the set sp ess H — U*=i [Aj nf i A^ up +4n] if n > 1. 

Let cosh -1 : [1, +oo) — > [0, +oo) refer to the function inverse to cosh : 
[0, +oo) — * [1, +oo), i.e., 



cosh 1 (i 



log( M + Vm 2 -1)- 



Further let TZ sl := \J r>1 W(K"). 



Theorem 18 Let w = e v be a weight in 7Z sl with MiRx—.oo v(x) = oo. Further 
let X be an eigenvalue of H such that A ^ sp ess H and assume that one of the 
following conditions is satisfied: 

(i) there is a j G {1, . . . , N} such that A G (A™ p + 4n, A"^) and 



lim sup 



dv(x) 



dx k 



< cosh 1 



min{A — X snp 



2n, Xfli - A + 2n} 



for every k = 1, . . . , n; 
(ii) X > A™ p + 4n and 



lim sup 



dv(x) 



for every k = 1, 
(ni) A < Ai nf an 



lim sup 



dx k 



dv(x) 



< cosh 1 



2n 



A - A s w up 



(24) 



2n 



2n 



, _i /Ai nf - A + 2n 

< cosh I — 

2n 



for every k — 1, . . . , n. 

Then every X-eigenf unction of H belongs to each of the spaces V '(IT 1 , C N , w) 
with p G (1, oo). 

Proof. For \i G [0, 1], let H' wli := w^H'w-^. The limit operators - XE 
are unitarily equivalent to the operator of multiplication by the matrix-function 



-2 Y cos(tpj + i/i0?) + 2n - X \ E + $ 9 



where 



t tp n ) G [0, 2tt]™ and 0* := lim 



dv{g{m)) 
dx 3 
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Note that 

9t(WJU(V0) = ^-2^cos^cosh/i^ ;i +2n- AJ E + $ 9 , 

where 6 9 W j := ^§^j ■ It is easy to check that condition (|24|) implies that 
A ^ specf/^ for every limit operator iJ^,,, of H w y. and every /j, € [0, 1]. Hence, 
by Theorem [T21 every A- eigenf unction belongs to l p (Z n ,C N ,w) for every 1 < 
p < oo. ■ 

Corollary 19 In each of the following cases 

{i) A S (A" up + An, Xf^) for some j G {1, . . . , N} and 

( min{A - \f P - 2n, - A + 2n}\ 
< r < cosh ^ — , 

V 2n J 

(ii) A > A™ p + An and < r < cosh" 1 ( A ~ A °£~ 2 " ) , 
(tit) A < Ai nf and < r < cosh" 1 ( A "' f ~* +2n ) , 

every \- eigenf unction of H belongs to ( p (Z",C iv ,e r ' 1 ") for each p G (l,co). 



Remark 20 in f/ie case o/ £/ie scalar Schrddinger operator (|£_?p wrf/i $ G 
SO{Z n ), we have 

sp ess = [$ inf , $ su p + An] 

with $ mf = liminf 2,^00 <&{x) and $ sup = limsupj,^^ // one o/ </ie follow- 

ing conditions holds for an eigenvalue A of H: 

(i) A > $ sup + 4n and < r < cosh" 1 ( X ~*'Z~ 2n ) » or 

(n) A < $ inf and < r < cosh" 1 f* int + 2 »-^ 



£/ien every X- eigenf unction of H belongs to l p (7j n , C N , e 1 *' 1 ') /or eachp £ (l,oo). 

4 The discrete Dime operator 

4.1 The essential spectrum 

On I 2 [2? , C 4 ), we consider the Dirac operators 

X> := Z> + e$I and V := cM fe 7 fe + c 2 m-f (25) 
where the 7 , k = 0, 1, 2, 3, refer to the 4x4 Dirac matrices, i.e., they satisfy 

TV + tV = 2o>£ 4 (26) 
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for all choices of j, k = 0, 1, 2, 3 where E4 stands for the 4x4 identity matrix. 
Further, 

d k :=I-V ek , k = 1,2,3 

are difference operators of the first order, h is Planck's constant, c the light 
speed, m and e are the mass and the charge of the electron, and <& is the electric 
potential. We suppose that the function $ is real-valued and belongs to the 
space SO {I?). 

It turns out that the operator V is not self-adjoint on Z 2 (Z 3 , C 4 ). Therefore 
we introduce self-adjoint Dirac operators as the matrix operators 

D := D + e&I with D : = ( ^° ) ' 

acting on the space Z 2 (Z 3 , C 8 ) (i.e. I refers now to the identity operator on that 
space). First we are going to determine the spectrum of D . It is 

(Bo-AI)(B +AI)= ( £( A) £( ° A) ) (27) 

where £(A) = h 2 c 2 Y + (m 2 c 4 - A 2 )/, and 

3 3 

k=l fc=l 

is the discrete Laplacian with symbol 

3 

f (f) = r(^i, <f2, <p 3 ) = ^2(2 - 2cos^ fe ), (p k £ [0, 2tt}. 
k=i 

Similarly, we denote by D (</?) and £(A, ip) the symbols of the operators Do an d 
C(X), respectively. Then 

(B (<p) - XE s )(B ((p) + XE 8 ) = C{X,<p)E s (28) 
with the scalar-valued function 

3 

l{X,ip) = h 2 c 2 J^(2 - 2 cos tp k ) + m 2 c A - A 2 . 

k=l 

We claim that A £ spec Do if and only if there exists a ipo € [0, 2tt] 3 such 
that £(A, ipo) — 0. Indeed, let A £ specD . Then there exists a tp £ [0,27r] 3 
such that det(D (^ ) - XE 8 ) = 0. Hence by (J2SJ C(X,(f ) = 0. Conversely, if 
£(A, <po) = 0, then it follows from (|28[) that 

(Doteo) - A£ 8 )(B (W)) + A£ 8 ) = 0. 
Hence, det(Do(<,oo) — XE$) = 0, whence A £ specDo. 
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Since the equation £(A, (f) = has two branches of solutions (spectral 
curves), namely 



\±(<p) = ±yh 2 c 2 T(ip) + m 2 c 4 , ip e [0, 2tt] 3 , 
the spectrum of Bo is the union 



spec B = [-Vl2ft 2 c 2 +m 2 c 4 , -mc 2 ] U [mc 2 , Vl2fr 2 c 2 + m 2 c 4 ]. 



Our next goal is to determine the essential spectrum of 1 



Do + e$7. All limit 



operators of B are of the form B 9 



e$ 9 / where $ 9 = lim,- 



is the partial limit of $ corresponding to the sequence g : N — > Z 3 tending to 
infinity. By what we have just seen, this gives 



spec J 



[e$ 9 - \J\2h 2 c 2 + m 2 c 4 , e$ 9 - mc 2 ] 



U [e$ 9 + mc 2 , e$ 9 + v / 12ft 2 c 2 + m 2 c 4 ] 
Since sp e;5S B = U g specB 9 we obtain 

sp ess B = [e$ inf - y/\2H 2 c 2 + m 2 c 4 , e$ sup - mc 2 ] 



inf 



U [e$ inf + mc 2 , e$ sup + y/l2h 2 c 2 + mV], 

where 

$ lnf := lim inf $(x) and <I> sup := limsup 3>(x). 

a;— »-oc a;— »oo 

In particular, if e($ sup - $ inf ) < 2mc 2 , then the interval (e$ sup - mc 2 , e<& 
mc 2 ) is a gap in the essential spectrum of B. 

4.2 Exponential estimates of eigenfunctions 

The following is the analog of Theorem [TBI 

Theorem 21 Let A ^ sp ess B 6e an eigenvalue o/B : /p(Z 3 , C 8 ) -> l p {Z 3 , C 8 ) 
lozfft p € (l,oo). Assume further that the weight w = e v is in lZ sl and that 
lim^^oo v(x) = oo. If one of the conditions 

(i) A e (e$ sup — mc 2 , e$ mf + mc 2 ) and, for every j — 1, 2, 3, 



lim sup 



dv(x) 



dxj 



< cosh 1 



2 c 4 - max{(e$ inf - A) 2 , (e$ sup - A) 2 } + 6h 2 



,(29) 



(ii) A > e$ sup + \/l2h 2 c 2 + m 2 c 4 and, for every j = 1, 2, 3, 



lim sup 



dv(x) 
dx-j 



< cosh 



_! ( (e$ sup - A) 2 - m 2 c 4 - 6ftV 



(30) 
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(Hi) A < e$ lnf — \/12h 2 c 2 + to 2 c 4 and, for every j — 1, 2, 3, 



lim sup 



dv(x) 



dxj 



1 /( e $i«f - A) 2 - m 2 c 4 - 6^ 2 c 2 



< cosh- ( * ) , (31) 



is satisfied, then every X-eigenfunction of the operator D belongs to the space 
Z P (Z 3 , C 8 , w) for each p S (1, oo). 

Proof. We will prove the assertion in case condition (i) is satisfied. The other 
cases follow similarly. Further, since the essential spectrum of D and the spectra 
of the associated limit operators do not depend on p, we can assume that p = 2 
in this proof. 

Let condition (|2U)) hold, and let A be an eigenvalue in the gap (e$ sup — 
mc 2 , e$ mf + mc 2 ) of the essential spectrum. In order to apply Theorem [TT1 to 
determine the decaying behavior of the associated eigenfunction u\, we need 
estimates of the spectrum of the limit operators (D w ^) 9 of D wf . :— it^Ow^' for 
/i S [0, 1]. The limit operator {w^Ve k w~^) 9 of w^V ek w~^ is of the form 

Hence, 

3 

(V W »Y = ^Tcj k {I - e^^i^T V ek ) + mc 2 j° + e<P g E 4 (32) 

k=\ 

where (&) ' = linw.- • 

Let W = D - e$I. The identity ((32]) implies that (D^ - AI)(B^ + A/) is 
the diagonal matrix diag (F, F) with 

F := ft 2 c 2 r^ + (m 2 c 4 - (e$ 9 - A) 2 )/ 

and 

=J2 [2I-e-(^Yv ek -ei^Y 

k=l ^ 

The operator is unitarily equivalent to the operator of multiplication by 
the function 

3 / / ( g v \ a 

v w»{v>) = f^(^i, ¥>2, ^3) = 51 2- 2cos iipk+i 



acting on the space L 2 ([0, 27r] 3 ). Note that 



k=l 

2/Tn 0^-13' 



<9.r 



3 / \ I 

9t(f» „ (^)) = 6 - 2 ^ cos ^ fc cosh ( j 



Hence, and by condition ([29]) . 

SK(n 2 c 2 f» „ fo>) + m 2 c 4 - (e$ 9 - A) 2 ) ^ (33) 
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for every sequence g denning a limit operator and for every fi 6 [0, 1]. The 
property ([33]) implies that A ^ specO^ M for every limit operator and every 
/x £ [0, 1]. By Theorem ITT1 every A- eigenf unction belongs to P(Z 3 , C 8 , vu) for 
every p £ (1, oo). ■ 

For the important case of the symmetric weight w(x) — e r ' x ', we obtain the 
following corollary of Theorem [^Tl 

Corollary 22 Let A be an eigenvalue of® : P{Z 3 , C 8 ) -> Z P (Z 3 , C 8 ). // one of 
the conditions 

(i) A € (e$ sup - mc 2 , e$ inf + mc 2 ) and 

_ x / m 2 c 4 + 6h 2 c 2 - max { (e$ inf - A) 2 , (e$ sup - A) 2 }\ 

°< r < C08h ^ ^2 J > 

(n) A > e$ sup + Vl2h 2 c 2 + m 2 c 4 and 

_! / (e$ sup - A) 2 - m 2 c 4 - 6^ 2 c 2 



< r < cosh 



6fi 2 c 2 



(Hi) A < e$ inf - \JVlh 2 c 2 + m 2 c 4 and 

_! / (e$ inf - A) 2 - m 2 c 4 - 6fi, 2 c 2 



< r < cosh 



V 6h 2 c 2 



is satisfied, then every X- eigenf unction of the operator ® belongs to the space 
P(Z 3 , C 8 , e^l) for every pe (1, oo). 

5 The square-root Klein-Gordon operator 

5.1 The essential spectrum 

Here we consider the square-root Klein-Gordon operator on Z 2 (Z n ), that is the 
operator 

K = V / c 2 h 2 T + m 2 c 4 + e$ 

where m > is the mass of the particle, fi. > is Planck's constant, c > the 
light speed, 4> e SO(Z n ) a scalar potential, and 

n 



is the discrete Laplacian on Z". The operator Kq :— \J c 2 h 2 T + m 2 c 4 is under- 
stood as the pseudodifference operator with symbol 



fc(r) = dc 2 H 2 T(T) + m 2 c 4 e S, 



19 



where f (r) = Y^=i( 2 - T j ~ T j 1 ) at r = (n, . . . , r„). Let 



f (?) := f(e^) = £(2 - 2cos^), y> = fa, • • • , € [0, 27r] n . 

3=1 

Every limit operator of K is unitarily equivalent to an operator of multiplication 
by a function of the form 



K B {v) = \J<?fr 2 T(<p) + m 2 c 4 + e$ 9 with $ s £ 



acting on L 2 ([0, 2tt]"). Thus, 



spec K 9 = (J [mc 2 + e$ 9 , V ' Anc 2 h 2 + m 2 c 4 + e$ 9 ] , 



where the union is taken with respect to all sequences g tending to infinity such 
that the partial limit $ 9 := lim„ woo $(g(m)) exists. Consequently, 



a Pess K = [mc 2 + e$ inf , V ' <±nc 2 h 2 + m 2 c 4 + e$ sup ]). 



5.2 Exponential estimates of eigenfunctions 

Theorem 23 Let A be an eigenvalue of the square-root Klein- Gordon operator 
K such that A ^ sp ess K, and let w = e° be a weight in lZ sl with lim^^oo v(x) = 
oo . If one of the conditions 



(i) A > e$ sup + \jAnh 2 c 2 + m 2 c A and 



lim sup 

X — >OG 

=,inf 



dv(x 



dx 3 



< cosh 



_! /m 2 c 4 - (e$ sup - A) 2 + 2n/iV 



(ii) A < e$ inf - yj4nh 2 c 2 + m 2 c A and 
dv(x) 



lim sup 



9a; i 



< cosh 



2 4 
_1 / TO C 



2nh 2 c 2 



(e$ inf - A) 2 + 2nh 2 & 



2nh 2 c 2 



(34) 



(35) 



is satisfied, then every X-eigenf unction of K belongs to Z P (Z", w) for every p £ 
(1, oo). 

Proof. The proof proceeds similarly to the proof of Theorem [2lJ It is based 
on the following construction. Let w = e v £ lZ sl . Then the limit operator 
is unitarily equivalent to the operator of multiplication by the function 



K^(ip) = ^Jc 2 h 2 T(ip + i(Vv))9 + m 2 c 4 + e$ 9 
acting on L 2 ([0, 2ir] n . Hence, 

£ 9 W „ (ip, A) := (Kl„ (<p) - A) (^]c 2 h 2 f{ip + i{Vv))9 + m 2 c i - (e$ 9 - A) 
= c 2 H 2 f(ip + i(\7v)) 9 + m 2 c 4 — (e$ s — A) 2 , 
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and 



(<p, A)) = c 2 H 2 ( 2 - costfj cosh(^)s ) + m 2 c 4 - (e$ 9 - A) 2 



Note that ^(Z^m ((/?, A)) ^ for every A satisfying condition (i) or (ii). Hence, 
A ^ sp ess K w n for every /it € [0, 1]. Thus, by Theorem [TT| every A-eigenfunction 
belongs to the space l p (Z n , w) for all p S (1, 00). ■ 

Specifying the weight in the previous theorem asw(i) = e r ' 1 ', we obtain the 
following. 

Theorem 24 Let A be an eigenvalue of K such that A ^ sp ess -ftT. // one 0/ the 
conditions 



(i) A > e$ sup + V4n^ 2 c 2 + to 2 c 4 ] and 

_! /m 2 c 4 - (e$ sup - A) 2 + 2nfr 2 c 2 



< r < cosh 



2nfi 2 c 2 



(«) A < e$ inf - \/Anh 2 c 2 + m 2 c 4 and 

_ a /m 2 c 4 - (e$ inf - A) 2 + 2nh 2 c 2 



< r < cosh 



2nh 2 c 2 



is satisfied, then every X-eigenf unction of K belongs to the space P(Z n , e r ' x ') 
for every p £ (1, 00). 
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